Graphene single electron transistor as a spin sensor for magnetic adsorbates 
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We study single electron transport through a graphene quantum dot with magnetic adsorbates. 
We focus on the relation between the spin order of the adsorbates and the linear conductance of the 
device. The electronic structure of the graphene dot with magnetic adsorbates is modeled through 
numerical diagonalization of a tight-binding model with an exchange potential. We consider sev- 
eral mechanisms by which the adsorbate magnetic state can influence transport in a single electron 
transistor: by tuning the addition energy, by changing the tunneling rate and, in the case of spin po- 
larized electrodes, through magnetoresistive effects. Whereas the first mechanism is always present, 
the others require that the electrode has either an energy or spin dependent density of states. We 
find that graphene dots are optimal systems to detect the spin state of a few magnetic centers. 

PACS numbers: 



I. INTRODUCTION 

Graphene is a very promising candidate for high pre- 
cision molecular sensing, due to its extremely large sur- 
face to volume ratio, and its electrically tunable large 
conductivity.— ~— On the other hand, being a zero-gap 
semiconductor with small mass and small density of spin- 
full nuclei, makes graphene a material with potentially 
large spin lifetime both, for carriers and host magnetic 
dopants^ Taken together, these two ideas naturally lead 
to the use of graphene as a detector of the spin state 
of extrinsic magnetic centers, in the form of magnetic 
adatoms, vacancies and spinfull molecules. This connects 
with recently reported^ experiments in which gated 
graphene nanoconstrictions, operating in the single elec- 
tron transport (SET) regime, showed hysteresis in the 
linear conductance when a magnetic field is ramped. This 
behavior was observed both when the molecular mag- 
nets were intentionally deposited on graphene, as well 
as in carbon nanotubes^ but also in the case of bare 
graphene nanojunctions^ where some type of graphene 
local moments^— is probably playing a role. 

The graphene spin sensor experiments of Ref. [H, are 
performed in the Coulomb Blockade regime, showing a 
vanishing linear conductance except in the neighborhood 
of specific values of the gate voltage V g . This means that 
the graphene nanoconstriction is weakly coupled to the 
electrodes and has a charging energy larger than the ther- 
mal energy (T ~ 100 mK). The height of the linear con- 
ductance peaks is significantly smaller than Gq = 2e 2 /h, 
the quantum of conductance. These conditions imply 
that transport takes place in the sequential regime*^ 
Thus, current flow takes place due to sequential tunnel- 
ing of electrons through the graphene constriction, which 
we refer to as the central region in the rest of the pa- 
per, and the entire device behaves like a single electron 
transistori n i 12 

The aim of this work is to provide a theoretical back- 
ground to understand how the magnetic state of local- 




FIG. 1: (Color online) (a) Scheme of a graphene constriction 
with randomly distributed magnetic centers, (b) Diagram 
with the system energy levels and graphene density of states. 



ized magnetic moments affects transport through the 
graphene nanoconstriction in the SET regime. This 
is different from previous works where the influence of 
the magnetic state of magnetic edgesi^ and adsorbed 
hydrogens^ on the conductivity was studied in the bal- 
listic regime, with a central island strongly coupled to 
the electrodes, and also in the diffusive regime,— as well 
as SET through graphene islands with magnetic zigzag 
cdgcsJ^ 
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The paper is organized as follows. In Sec. [TT| we dis- 
cuss a tight-binding Hamiltonian for the graphene island 
exchanged coupled to the spins of the magnetic adsor- 
bates. The results of this microscopic calculation justify 
the use of a simple single-orbital spin-split model for the 
SET, discussed in section HTT1 together with the possible 
mechanisms that enable the magnetic sensing. Finally, 
conclusions are presented in Sec. IIVI 

II. MODEL FOR GRAPHENE ISLAND WITH 
MAGNETIC ADSORBATES 

A. Hamiltonian 

Our starting point is a microscopic model for elec- 
trons confined in a graphene nanoisland which are ex- 
changed coupled to the magnetic centers. The graphene 
central island is described with a tight-binding Hamilto- 
nian for the honeycomb lattice contained in a rectangular 
stripe of dimensions L x x L y . In order to avoid the spin- 
polarized states formed at the the zigzag edges) 16 ' 17 we 
impose periodic boundary conditions in one direction so 
that the structure only has open edges of armchair type. 
The coupling to the magnetic moments of the adsorbate 
molecules m(i) is then assumed to be a local exchange J 
or spin-dependent potential, affecting N sites randomly 
selected in the graphene central island. For simplicity, 
we consider that the magnetic moments of the molecules 
are all oriented along the same axis, which we choose as 
the spin quantization axis. These assumptions are good 
approximations in the case of strongly uniaxial TbPC2 
molecules^ Hence, we can write the Hamiltonian of the 
graphene and adsorbates as 

U = n a + J$^ro,(i)S,(i) + eV g (n tot -ft), (1) 

i 

where Ho is the tight-binding Hamiltonian for 7r-electrons 
in graphene considering nearest neighbor interactions, J 
is the strength of the exchange coupling between the 
graphene electrons and the magnetic moment of the 
molecules, which can take two values, m z (i) = ±1. Fi- 
nally, the last term in the Hamiltonian describes the elec- 
trostatic coupling of the total charge of the dot, which 
can be either or e, given by the difference in the number 
of electrons N and the number of carbon atoms iVxoT in 
the central island. S z (i) is the local spin density of the 
p z electrons in graphene at site i, 

S z (i) = - (4 t Cit - ct^) , (2) 

where creates an f electron at the p z orbital of site i 
of graphene. In the following we assume that magnetic 
fields controlling the spin orientation of the adsorbates, 
are applied along the plane of graphene so that it is a 
good approximation to neglect the diamagnctic coupling 
to the graphene electrons. 



There are several independent microscopic mechanisms 
for spin dependent interaction between magnetic adsor- 
bates and the graphene 7r electrons that can be modeled 
with Eq. ([T|). In the case of magnetic molecule such as 
TbPC2, used in in Ref. H, 

the magnetic Tb atom is sep- 
arated from the graphene electrons by the non-magnetic 
atoms of the molecule, and the most likely mechanism 
for spin coupling is kinetic exchanged This coupling 
will generate a local Kondo-exchange between graphene 
electrons and the molecules More complicated scenar- 
ios, like coupling of graphene electrons to unpaired elec- 
trons in the organic rings of the molecules, would imply 
that every molecule affects several sites in graphene. Di- 
rect dipolar coupling would also affect several sites per 
molecule, but the average magnetic field created by a 
magnetic moment of 5 \xb at 0.5 nm on a disk with an 
area around 400 nm 2 , the graphene constriction area in 
Ref. [U, is smaller than 1 //T, which would yield a negli- 
gible maximal Zeeman coupling of neV per molecule. 

B. Relevant energy scales 

The reported dimensions of the central region, L x ~ 
L y ~ 20 nm, lead to an energy spacing Ae of the sin- 
gle particle spectrum much larger than the temperature 
and the charging energy n i 12 We also assume that the ex- 
change induced shifts are smaller than the single particle 
splitting. As a result, the effect of exchange is to shift 
the bare energy levels, without mixing them. Thereby, 
we can safely assume that electrons tunnel through just 
one of the single particle levels, which might be spin- 
split due to exchange with molecules. We assume that 
the charge of the central island fluctuates between q = 
and q = — \e\, and that the transport level is the low- 
est unoccupied level of the central island spectrum. The 
energy of the transport level reads: 

^T = e +^-\e\V g , (3) 

where eo is the single particle electron level, a = ±1 de- 
notes the spin direction and S is the magnitude of the spin 
splitting, which is a functional of the magnetic landscape 
{m z }. 

Within our model, a given magnetic landscape is de- 
fined by the location i and the magnetic state m z (i) of 
the N magnetic adsorbates. In Figs. [2] (a, b), we plot the 
value of 5 for all the possible magnetic states of a given 
arrangement of N = 10 atoms, for two different values 
of J. This choice corresponds to the estimated number 
of molecules in Ref. 0. These figures show a correlation 
between the magnitude of the splitting S and the total 
magnetization Mr. The dispersion of 6 for a fixed total 
magnetization Mr is the outcome of indirect exchange 
coupling^ 

For comparison with the experiments, it is worth con- 
sidering two extreme magnetic landscapes. At large ex- 
ternal field, all the magnetic moments are aligned, i.e., 
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FIG. 2: Spin splitting S for a sample of dimensions L y = 
N y a and = N x y/2>a, with AT^ = 15 and N y = 17, being 
a = 2.46A the lattice constant of graphene. 8 versus total 
magnetization Mt for a single spatial distribution of N = 10 
magnetic impurities with J = 5 meV (a) and J = 10 (b). 
The red solid line indicates the average value {S) p0 s- Average 
(5)™ versus number of magnetic molecules N for J = 5 meV 
(c) and exchange energy J for A^ = 4 (d). The error bars 
correspond to the standard deviation. 



m z (i) = +1. We refer to this as the ferromagnetic (FM) 
landscape. At magnetic fields smaller than the coercive 
field of the magnetic molecules, their average magnetiza- 
tion should be zero and thus, X)i m z(*) = 0- We refer 
to these cases as non- magnetic (NM). In order to sample 
the positional disorder we perform an average over posi- 
tional configurations, both for NM and FM cases. For a 
fixed spin choice {m z } with Mt = 0, an average over po- 
sitional configurations yields (<5) pos = 0. The reciprocal 
statement is also true: for a fixed positional configura- 
tion, an average over all the magnetic landscapes with 
Mt = also yields an average (<5) sp ms = 0. 

In Fig. |we plot the average <<5}™ - (S)™ = (S)™ 
over 500 realizations function of the number of 

molecules N [Fig. EJc)] and as a function of the molecule 
electron exchange J [Fig. E^d)]. We have also calcu- 
lated (<5)™ fixing the number of magnetic centers TV, the 
strength of the coupling J and changing A^tot = N x N y , 
the total number of carbon atoms in the island. We find 
that the results of all these simulations can be summa- 
rized in the following equation: 



(S) 



FM 
pos 



N 



Ntot 



(4) 



Whereas this result has been obtained from exact numer- 
ical diagonalization of the Hamiltonian, this dependence 
can be rationalized using first order perturbation theory, 
which yields the spin dependent shift of the transport 
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(5) 



where <j>T{i) is the J — wave function of the transport 
orbital. We now use |0t(«)| 2 — Nt t so ^ a t we can 
approximate: 
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(6) 



Using the fact that Mt = N for the FM configurations 
and for the NM ones, we arrive to Eq. (j4|). 



III. SPIN-SPLIT SINGLE ORBITAL MODEL 
FOR SET 

In this section we discuss SET across a central island 
with a single spin split particle level. This is justified 
by the results of the previous section. We obtain expres- 
sions for the current of the system and we discuss the 
conditions under which the conductance depends on the 
magnetic state of the single electron transport. 



A. Single electron transistor with a spin-split 
single orbital model 

We consider single electron transport though a spin 
split single transport level^ with energy e T - We as- 
sume that the occupation of the transport level can be 
either or 1, the doubly occupied configuration being 
much higher in energy. Within these approximations, 
the transport level has three relevant many body states: 
uncharged, and the two charged with ^ or J, spins. In 
the zero-applied bias limit, each of these states will be 
occupied according to the thermal equilibrium distribu- 
tion, which we denoted as Pq, P-\ and Pj, respectively. 
In the SET regime we are interested and within the lin- 
ear response (eVuas *C ksT), transport will be enabled 
only when the addition energy lies within the thermally 
broadened transport window defined by the applied bias. 

Under these approximations, the current flowing from 
the left electrode to the central island is given by: 



(7) 



where W^ CT and W^ are rates for electron tunneling 
from the left electrode to the dot and vice-versa. Conti- 
nuity equation ensures that this current is identical to the 
current flowing towards the right electrode and, thereby, 
equal to the net current flow. The tunneling rates for 
electron tunneling out of and into the dot^ are given 
respectively by 
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- -jr \Tl\ pl(^t) [1 - /(A. + Vuas)] , (8) 
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and 

Wtf-a = t |T l | 2 P l(4)/(A. + H ios ), (9) 

where Tl is the strength of the dot - left electrode cou- 
pling, and 

A a EE e£ - £ F , (10) 

are the spin dependent addition energies. Importantly, 
both S and appear on equal footing, as additive quanti- 
ties in this equation. The density of states of the left elec- 
trode, evaluated at the spin-dependent transport level en- 
ergy, is denoted by pi(ej.), while f(e a ) = (e^ A " + l) 
denotes the Fermi function. The electrode Fermi energy 
Ep is taken to change linearly with the bias voltage Vbias- 
In the zero bias limit, the linear conductance reads: 

G = Go £ (P + P a ) ^Sech 2 , (11) 

where 

T a /h=^\T L f PL (e° T ), (12) 

is the single particle tunneling rate between the electrode 
and the transport level. 

B. Influence of magnetic state on conductance 

From the above discussion it is apparent that, for a 
given gate potential V g and temperature T, the linear 
conductance depends on the magnetic landscape affect- 
ing the central island through two classes of independent 
mechanisms, illustrated in Fig. [3j 

1. The change of the addition energies A CT which, as 
we show below, would result in a lateral shift of the 
G(V g ) resonance curve [Fig. Ufa)]. 

2. The change of the electron lifetime T = J2a^<?' 
that would result in a vertical resizing of the G(V g ) 
resonance curve [Fig. GDJa)]. 

In the first mechanism the change in the magnetic 
state modifies the value of 5, which must have a simi- 
lar effect than changing the gate potential. It resembles 
the magneto- Coulomb effect ) 23 i 24 by which the applied 
magnetic field changes the Fermi energy of the electrode, 
shifting the G(V g ) curves. However, this first mechanism 
necessarily implies a change of sign of the variation of G 
as the gate potential is scanned along the resonance (see 
top panel of Fig. SJ. Importantly, this is not observed 
in the experiments with magnetic molecules^ but it is 
observed in the case of graphene nanoconstrictions.— 

Motivated by the behavior reported in Ref. [U, we pay 
attention also to the second mechanism. For spin un- 
polarizcd transport, the change in the transport energy 




FIG. 3: (Color online) (a) Scheme showing the transport level 
energy splitting. Scheme of spin-dependence of transport due 
to: (b) detuning of the transport level with respect to the 
electrode Fermi energy, (c) magnetoresistance associated to 
spin polarized electrode(s), and (d) energy dependent tunnel- 
ing rates. 

level results in a change on tunneling rate T/H only if the 
electrode density of states depends on energy, which is 
exactly the case of graphene. For spin-polarized trans- 
port, the relative orientation of the electrode and island 
magnetic moment give rise to magneto-resistive effects 
that are accounted for by the changes in F^ 



C. Transport for constant tunneling rates 

We now discuss our transport simulations for the 
graphene single electron transistor spin sensor. We focus 
on the first spin sensing mechanism in a single electron 
transistor: changes in spin splitting of the transport level 
produce changes in addition energies A CT (figure Qb). 
For that matter, we neglect both the energy and spin 
dependence of the tunneling rates T a . In Fig. Ufa) we 
show the linear conductance, in units of go = @TGq, as a 
function of gate voltage, for several values of the trans- 
port level splitting S, in units of fc^T. It is apparent 
that the Coulomb Blockade peaks undergoes a lateral 
shift, as expected from the fact that V g and S appear on 
equal footing on the spin dependent addition energies. 
At A = the two spin channels contribute. Therefore, 
as we increase 5, the height of the conductance peaks de- 
creases, because one of the two spin channels is removed 
from the transport window of width fc^T 

In figure (Qb we plot the variations in the linear con- 
ductances as a function of the spin splitting S, for several 
values of V g . We see two types of curves. For values of V g 
such that the transport level is occupied, as we increase 
(5 the transport level is pushed downwards, away from 
the elastic transport window, switching off the transistor 
conductance. In contrast, for V g such that the transport 
level is empty for 6 = 0, lying above the elastic transport 
window, ramping 6 makes one of the two spin states of 
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FIG. 4: (Color online) Conductance in units of go = Gof3T 
as a function of the gate voltage eV g and the level splitting 
S. Panel (a) shows the normalized conductance as a function 
of the gate potential e/3V g , where the labels corresponds to 
the different j35 values. Panel (b) shows the conductance as 
a function of /3S for several values of the gate e/3V ff . For the 
sake of clarity all curves have been displaced by 0.1. In (c, 
d) we present a contour plot of the dot charge (defined as 
Q — P-f + Pj,) and the dot magnetization (m = Pf — Pj_) as a 
function of the gate voltage eV g and level splitting 5. 



the transport level enter the transport window, giving 
rise to the double peak structure. The fact that S and V g 
play analogous roles is illustrated in Figs.[4^c, d), where 
we show the average magnetization and occupation of the 
transport level in the phase diagram defined by these two 
variables. 



D. Transport with energy dependent tunneling 
rates 

We now consider the second mechanism for spin sens- 
ing in a single electron transistor: changes in spin split- 
ting of the transport level produce changes in the tun- 
neling rates r,j. This can happen for two reasons: 

1. One of the electrodes is spin polarized, so that 
TV 7^ IV Spin polarized transport is sensitive to 
the product of the magnetic moment of electrode 
and central island. This type of effect has been 
thoroughly discussed in the case of SET with fer- 
romagnetic electrodes ! 25 ' 26 

2. The density of states of the electrode depends on 
energy. Thus, changes in the transport level change 
To-, for both spins. This is a natural scenario for 
graphene electrodesj 32 ' 33 

Let us consider first the case of spin polarized elec- 
trodes. We do the assumption that the density of states 
are spin dependent but energy independent: p a = po(l + 
aV), with V the electrode polarization. In Fig. [5ja) we 
plot the linear conductance vs V g curves for several values 
of the splitting S, assuming a large value of the electrode 
spin polarization V = 0.9. It is apparent that, on top 
of the shift of the resonance curve whose origin was dis- 
cussed in the previous subsection, there is a change in 
the amplitude of the curve. Notice that G(S = fcgT) 
and G(5 = 2ksT) are smaller than G(0) for all values 
of V g . In this specific sense, the gate- independent spin 
contrast is similar to the experimental report with mag- 
netic molecules^ In Fig. [Ha) we show the linear conduc- 
tance as a function of S for different values of V g . It is 
apparent that, as opposed to the case of non-magnetic 
electrodes shown in Fig. 5(b), the function G(5) is no 
longer an even functions, reflecting the magneto-resistive 
behaviour. Basically, transport is favored when the spin 
polarization of the electrode and the central island are 
parallel. 

We now consider a non-spin polarized electrode with 
an energy dependent density of states. This scenario 
occurs naturally in graphene. If we consider idealized 
graphene electrodes, neglecting effects of interactions, 
disorder and confinement, we have p(e) = po|e — Eo\, 
where Ed is the Dirac point. The G(V g ) curves, shown 
in Fig. [Hb) for different values of S, shift and change 
amplitude. The shift is related to the change of the ad- 
dition energies, discussed in the previous subsection, and 
the change in amplitude comes from the variation of the 
tunneling rate as the transport level scans the energy 
dependent density of states of the electrode. 

In Fig. [Hd) we plot G(6) for several values of V g . The 
curves are similar to the case with energy independent 
tunneling rates, except for the dip at zero S which occurs 
because we chose the bare transport level right at Dirac 
point. This is the most favorable choice to maximize the 
effect of energy dependence of T. From our results, and 
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FIG. 5: (Color online) Normalized Conductance for ferromag- 
netic electrodes as a function of the level splitting S for several 
gate values. In (a, c) the density of states depends of the po- 
larization p = po(l + <jV), in this particular case we take a 
polarization of V = 0.9. In (b, d) the electrode density of 
states is linear with the energy, p(A CT ) = p(eo + u5/2 — Ed)- 
For the sake of clarity the conductance curves in (c, d) have 
been displaced by 0.2 and 4 units respectively. 

given the fact that experimentally is not possible to put 
the Fermi energy arbitrarily close to the Diract point,— 
we find it unlikely that this effect is playing a role in the 
experiments. 

E. Sensitivity of the single electron transistor spin 
sensor 

We now discuss the sensitivity of the spin sensor based 
on the graphene single electron transistor, as described 
by our model, neglecting changes in T. From Fig. (Ha) we 
propose, as rule of thumb, that variations of 8 similar or 
larger than ksT can be resolved. Estimating 8 from the 
case of fully spin polarized magnetic adsorbatcs, given in 
Eq. @, we find a relation between the minimal num- 
ber of magnetic centers N that can be detected, and the 
temperature and exchange constant: 



It is apparent that decreasing the temperature, or in- 
creasing the spin-graphene exchange coupling, increases 
the sensitivity of the device (makes it possible to detect 
a smaller concentration of molecules). For instance, at 
100 mK, and taking -/Vtot ~ 15000, which corresponds 



to an approximate area of 400 nm 2 , one could detect 10 
molecules for an exchange coupling J > 15 meV. 

Recent reports have shown that it is possible to fab- 
ricate graphene nano islands with lateral dimensions of 
lnm.— These dots have -/Vtot < 100. Thus, they would 
permit the detection of the spin of a single magnetic ad- 
sorbate provided that ksT is kept hundred times smaller 
than J. For T = 100 mK this implies J > 1 meV. Inter- 
estingly, in such a small dot Coulomb Blockade persists 
even at room temperature, but increasing ksT keeping 
the sensitivity would require also to increase J. 



IV. DISCUSSION AND CONCLUSIONS 

Because of its structural and electronic properties, 
graphene is optimal for a spin sensor device. Being all 
surface, the influence of adsorbates on transport should 
be larger than any other bulk material. Because of the 
linear relation momentum and large Fermi velocity, en- 
ergy level spacing in graphene nano structures can easily 
be larger than the temperature, the tunneling induced 
broadening, and the perturbations created by the adsor- 
bates. One of the consequences is that single electron 
transport takes place through a single orbital level. 

Our simulations show how the spin splitting 8 of the 
transport level is sensitive to the average magnetization 
of the magnetic adsorbates, which is controlled by appli- 
cation of a magnetic field along the plane of graphene, to 
avoid diamagnctic shifts. On the other hand, the linear 
conductance G of the single electron transistor depends 
on (5, which accounts for the sensing mechanism. More 
specifically, G depends on 8 due to either changes in the 
spin dependent addition energies A CT or changes in the 
electrons lifetime T a . The first is independent of the na- 
ture of the electrodes, whereas the second only happens 
if they are magnetic or have an energy dependent density 
of states. 

We have shown how, within an independent particle 
model and in the single electron transport regime, the 
energy dependence of the graphene electrode density of 
states can only be relevant if the transport energy level is 
fine tuned to the Dirac point. However, this fine tuning 
is quite unlike to happen in experimental conditions 2 ^. 
Still, the combined action of disorder and Coulomb in- 
teraction could give rise to a so called Coulomb Gap in 
the density of states of graphene, that might make the 
tunneling rates depend on the energy^ - — 

Finally, we have assumed that both the edges of 
the graphene island and graphene electrodes are non- 
magnetic. Our discussion of the effect of spin-polarized 
electrodes on the transport properties of the device would 
be valid for electrodes with ferromagnetic zigzag edges. 
A second possibility, out of the scope of this work, is to 
consider a graphene single electron transistor whose cen- 
tral island has ferromagnetic edges. This case has been 
already studied.— 

In conclusion, we have studied the mechanisms by 
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which a graphcnc single electron transistor could work 
as a sensor of the magnetic order of magnetic atoms or 
molecules adsorbed on the graphene central region. Our 
work has been motivated in part by recent experimen- 
tal works'. Whereas further work is still necessary to 
nail down the physical mechanisms for the spin sensing 
principles underlying the experimental work, our study 



provides a conceptual framework for graphene single elec- 
tron transistor spin sensors. 
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